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(54) Method for determining the road handling of a tyre of a wheel for a vehicle 



(57) A method for determining the road handling of 
a tyre, comprising descriptions of the tyre by means of 
a first, concentrated-parameter, physical model (1 ) and 
by means of a second, finite-element mode! (30). a sim- 
ulation on the second, finite-element model (30) of a se- 
lected series of dynamic tests and an application to the 
first physical model (1 ) of equations of motion suitable 
for representing the dynamic tests in order to obtain first 
and second frequency responses of selected quantities; 
a comparison between the first and second frequency 
responses of the selected quantities for determining the 
concentrated parameters of the first physical model (1) 
and physical quantities indicative of the drift behaviour 
of the tyre (Fig. 1 ). 
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Description 

[0001] This invention relates to a method for determining the road handling of a tyre of a wheel for a vehicle. 
[0002] At the present time, to determine the road handling performance of a tyre, the manufacturers of pneumatic 
s tyres are obliged to produce numerous physical prototypes in order to experimentally evaluate the effects of the various 
design parameters on the drift behaviour of the tyre, under steady state and transient state conditions. The experimental 
tests are conducted according to iterative procedures, that are largely empirical and based on experience and are also 
extremely demanding in terms of time and cost. 

[0003] Furthermore, automobile manufacturing companies are insisting more and more frequently that the makers 
10 of pneumatic tyres come up with tyres with extremely precise technical characteristics as early as the initial stages of 
vehicle study and dynamic behaviour forecasting. 

[0004] In such a position, the tyre manufacturers are finding it very difficult to respond satisfactorily and with the 
necessary flexibility to the various market demands. 

[0005] The object of this invention is to provide a scientific methodology with which to identify the performance char- 
ts acteristics of a tyre in relation to road handling, on the basis of previously defined design specifications. 

[0006] The above object is achieved according to this invention by a method for determining the road handling of a 
tyre of a wheel lor a vehicle, said tyre being made from selected mixes of rubber and reinforcing materials, said method 
comprising: 

20 a) a first description of said lyre by means of a first, concentrated-parameter, physical model, said first physical 

model compiising a ngid ring which represents the tread band provided with inserts, a belting structure and cor- 
responding carcass portion of said tyre, a disk which represents a hub of said wheel and beading of said tyre, 
principal springs and dampers connecting said rigid ring to said hub and representing sidewalls of said tyre and 
air under pressure inside said tyre, supplementary springs and dampers representing deformation phenomena of 

25 said belting structure through the effect of a specified vertical load, and a brush model simulating physical phe- 

nomena in an area of contact between said tyre and a road, said area of contact having a dynamic length 2a, 

b) a definition of selected degrees of freedom of said first physical model, and 

c) an identification of equations of motion suitable for describing the motion of said first physical model under 
selected dynamic conditions, 

30 characterized in that it comprises 

d) the definition of said concentrated parameters, said concentrated parameters consisting of the mass M c and a 
diametral moment of inertia J c of said rigid ring, the mass M m and a diametral moment of inertia J m of said disk, 
structural stiffnesses K c and structural dampings R c respectively of said principal springs and dampers, and residual 
stiffnesses K r and residual dampings R r respectively of said supplementary springs and dampers, wherein 

35 

said structural stiffnesses K c consist of lateral stiffness K cy between said hub and said belt, camber torsional 
stiffness K c6x between said hub and said belt and yawing torsional stiffness K cez between said hub and said belt, 
said structural dampings R c consist of lateral damping between said hub and said belt, camber torsional 
damping R c9x between said hub and said belt and yawing torsional damping R c9z between said hub and said 
40 belt, 

said residual stiffnesses K r consist of residual lateral stiffness K^, residual camber torsional stiffness and 
residual yawing torsional stiffness K^, and 
' - said residual dampings R r consist of residual lateral damping ft^, residual camber torsional damping R^ and 
residual yawing torsional damping R^, 

45 

e) a description of said tyre by means of a second, finite-element model comprising first elements with a selected 
number of nodes, suitable for describing said mixes, and second elements suitable for describing said reinforcing 
materials, each first finite element being associated with a first stiffness matrix which is determined by means of 
a selected characterization of said mixes and each second element being associated with a second supplementary 

so stiffness matrix which is determined by means of a selected characterization of said reinforcing materials, 

f) a simulation on said second, finite-element model of a selected series of virtual dynamic tests for exciting said 
second model according to selected procedures and obtaining transfer functions and first frequency responses of 
selected quantities, measured at selected points of said second model, 

g) a description of the behaviour of said first physical model by means of equations of motion suitable for repre- 
ss senting the above dynamic tests for obtaining second frequency responses of said selected quantities, measured 

at selected points of said first physical model, 

h) a comparison between said first and said second frequency responses of said selected quantities to determine 
errors that are a function of said concentrated parameters of said first physical model, and 
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i) the identification of values for said concentrated parameters that minimize said errors so that said concentrated 
parameters describe the dynamic behaviour of said tyre, 

j) the determination of selected physical quantities suitable for indicating the drift behaviour of said tyre, and 
k) the evaluation of the drift behaviour of said tyre by means of said physical quantities. 

5 

[0007] To advantage, said selected physical quantities are the total drift stiffness K d of said tyre, in turn comprising 
the structural stiffness K c and the tread stiffness K b , and the total camber stiffness Ky of said tyre. 
[0008] According to a preferred embodiment, said method also comprises 

io |) a definition of said brush model, said brush model having a stiffness per unit of length c py and comprising at 

least one rigid plate, at least one deformable beam having a length equal to the length 2a of said area of contact 
and at least.one microinsert associated with said beam, said microinsert consisting of at least one set of springs 
distributed over the entire length of said beam, said springs reproducing the uniformly distributed, lateral and 
torsional stiffness of said area of contact. 

15 

[0009] Preferably, said degrees of freedom referred to at previous point b) are composed of: 

absolute lateral displacement y m of said hub, absolute yaw rotation a m of said hub and~absolute rolling rotation 
P m of said hub, 

20 - relative lateral displacement y c of said belt with respect to said hub, relative yaw rotation a c of said belt with respect 
to said hub and relative roiling rotation p c of said belt with respect to said hub, 

absolute lateral displacement y b of said plate, absolute yaw rotation o b of said plate and absolute rolling rotation 

P b of said plate, 

and 

25 - absolute lateral displacement y s of the bottom ends of said at least one microinsert. 

[0010] According to another embodiment, said selected series of virtual dynamic tests referred to at previous point 
f) comprises a first and a second test with said tyre blown up and not pressed to the ground, said first test consisting 
in imposing a translation in the transverse direction y on the hub and in measuring the lateral displacement y c of at 

30 least one selected cardinal point of said belt and the force created between said hub and said belt in order to identify 
said mass M c said lateral stiffness K^, and said lateral damping R^, said second test consisting in imposing a camber 
rotation B x on said hub and in measuring the lateral displacement of at least one selected cardinal point of said belt y c 
and the torque transmitted between said hub and said belt in order to identify said diametral moment of inertia J c , said 
camber torsional stiffness K^, said camber torsional damping R c9x , said yawing torsional stiffness K cez and said 

35 yawing torsional damping R^. 

[001 1] Preferably said selected series of virtual dynamic tests referred to at previous point f) also comprises a third 
and a fourth test with said tyre blown up, pressed to the ground and bereft of said tread at least in said area of contact, 
said third test consisting in applying to said hub a sideward force in the transverse direction F y and in measuring the 
lateral displacement y c of said hub and of at least two selected cardinal points of said belt in order to identify said 

40 residual lateral stiffness K^, said residual lateral damping R^, said camber residual stiffness K^ XJ and said camber 
residual damping R rt)x , said fourth test consisting in applying to said hub a yawing torque C ez and in measuring the 
yaw rotation of said hub a m and the lateral displacement y c of at least one selected cardinal point of said belt in order 
to identify said residual yawing stiffness and said residual yawing damping R^ z . 
[0012] According to another embodiment, said method also comprises 

45 

m)an application to said first physical model of a drift angle a, starting from a condition in which said at least one 
beam is in a non-deformed configuration and said brush model has a null snaking o b , 

n) the determination of the sideward force and the self-aligning torque that act on said hub through the effect of 
said drift a and which depend on the difference a - o b and on the deformation of said at least one beam, 
so o) the determination of the deformation curve of said at least one beam, 

p) an application of said sideward lorce and said self -aligning torque to said second, finite-element model in order 
to obtain a pressure distribution on said area of contact and 

q) the determination of the sideward force and the self-aligning torque that act on said hub through the effect of 
said drift on said first physical model, that depend on the pressure distribution calculated in the previous step p), 
55 r) a check, by means of said pressure distribution obtained in the previous step p), that said sideward force and 

said self-aligning torque are substantially similar to those calculated in previous step q), 
s) a determination of the sideward force and of the self-aligning torque for said angle of drift, and 
t) repetition of the procedure from step m) to step s) for different values of the drift angle a to obtain drift, force 
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and self-alignment torque curves, suitable for indicating the drift behaviour under steady state conditions of said 
tyre, and 

u) the evaluation of the steady state drift behaviour of said tyre. 
5 [0013] According to another preferred embodiment, said method also comprises: 

i) a simulation of the behaviour of said first physical model in the drift transient state by means of equations of 
motion reproducing selected experimental drift tests, and 

ii) the determination,, with a selected input of a steering angle imposed on said hub, of the pattern with time of the 
io selected free degrees of freedom of said first physical model, of the sideward force and of the self -aligning torque 

in said area of contact in order to determine the length of relaxation of said tyre. 

[0014] To advantage, said first elements of said second, finite-element model have linear form functions and their 
stiffness matrix is determined by means of selected static and dynamic tests conducted on specimens of said mixes, 
is whereas the stiffness matrix of said second elements is determined by means of selected static tests on specimens 
of said reinforcing materials. 

[0015] With the method according to this invention, three main results are obtained: 

1 . determination of the links between physical parameters of the tyre and its structural properties; 
?0 2. determination of the steady state drift behaviour of the tyre, without the need to build prototypes at this stage; 

3. determination of the transient state behaviour of the tyre, when a generic law of motion is imposed on the hub, 
without the need to build prototypes at this stage. 



[0016] These results have been achieved by the production of a very simple, first physical model, with only nine 
25 degrees of freedom, that manages to make allowance for the majority of the structural characteristics of the actual tyre. 
[001 7] The structural characteristics of the tyre are reproduced in the first physical model by means of an appropriate 
condensation of concentrated equivalent masses, stiffnesses and dampings. 

[001 8] In practice, the concentrated-parameter model is equivalent to a kind of dynamic concentration of the complex 
finite-element model, summarizing all its dynamic characteristics in a low number of mass, damping and stiffness 
30 parameters. 

[0019] More specifically, it has been proven that this correspondence may be held valid in the range of frequencies 
between 0 and 80 Hz. 

[0020] The method enables identification of the structural parameters needed for the complete description of the first 
physical model using simulated virtual numerical tests with a second, extremely detailed model built from finite element 
35 models (F.E.M.) reproducing the behaviour of the non-rolling tyre (not drifting). 

[0021] One of the major advantages of the method according to the invention is that it partly dispenses with the need 
to construct physical prototypes and the resultant experimental tests in the iterative process of tyre determination, 
replacing this approach with virtual prototyping. 

[0022] The tyre's design parameters (characteristics of the mixes, inclination of the threads, shape of the sidewalls, 
40 width of the belt, etc.) are directly fed into the second, finite-element model which is extremely detailed. 

[0023] The concentrated parameters of the first physical model are identified by minimizing the difference between 
the vibrational dynamic behaviour of the second, finite-element model of the non-rolling tyre and the corresponding 
response given by the first physical model. 

[0024] The identification procedure defined comprises various operations that are executed in a precise, pre-estab- 
45 lished order. Starting from the transfer functions obtained by means of a series of virtual dynamic tests conducted on 
the non-rolling finite-element model, the masses, stiffnesses and dampings of the concentrated-parameter model are 
determined, providing a better description of the dynamic behaviour of the tyre. 

[0025] Thus the identification procedure enables a link to be established between the design parameters (fed into 
the second finite-element model) and the condensed structural properties (contained in the first model with nine degrees 

50 of freedom), which is extremely useful in the construction of a tyre. 

[0026] The method consists in linking the design parameters of the tyre, characteristics such as the mixture and the 
belt, to the structural parameters, for instance the structural stiffness and camber stiffness of the tyre, because the 
quantities appearing in the model used have a physical significance. This means that these quantities are directly 
linked to the design parameters, in other words the model used is a physical model. By so doing, any change to the 

55 design parameters of the tyre leads to a change of the parameters of the predictive physical model of the tyre and this 
change, in turn, produces a variation of the tyre's structural parameters. 

[0027] The model permits identification of the structural parameters starting from dynamic analysis made on the 
second, finite-element model of the non-rolling tyre. One requirement of the concentrated-parameter model is, in fact, 
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that it be predictive of actual behaviour of the tyre. 

[0028] One of the main advantages of the method according to the invention is that the concentrated parameters 
are not identified by means of experimental tests on prototypes, but by means of virtual. dynamic tests on the' finite- 
element model of the non-rolling tyre. 
5 [0029] The method according to the invention uses a model of contact between tyre and road that enables forecasting 
of the drift curves at steady state. Also fed into the brush model, in addition to the longitudinal and transversal stiffness 
of the inserts of the tread, was their torsional stiffness, these stiffnesses being identified by means of numerical sim- 
ulations on the second finite-element model, without any need for experimental tests. 
[0030] The method .according to the invention enables drift curves to be determined 

10 

by applying a drift angle a to the first physical model; 
- created in the area of contact on account of the drift are a sideward force and a self-aligning torque that act on the 
first physical model as forces acting on the hub and cause a lateral displacement and a snaking motion of the plate 
of the brush model; 

is - because of the snaking motion a b and the lateral displacement of the plate, the forces set up in the contact area 
are modified; this results in a variation of the unconstrained degrees of freedom, among which those of the plate, 
and therefore of the forces acting of the first physical model. 

[0031] With this procedure, after a certain number of iterations, a point is reached at which the degrees of freedom 
20 of the first physical model settle about a steady state value. In this situation, the sideward force and the self-aligning 
" torque created in the area of contact and from which the drift curves may be obtained are determined. 
[0032] The method according to the invention also enables transient state drift behaviour of the tyre to be evaluated, 
making allowance in the brush model for the dynamic deformations undergone by the inserts of the tread in this stage. 
[0033] In this way, the length of relaxation of the tyre while drifting is determined upon variation of the running con- 
25 ditions (speed, vertical load, drift angle, etc.). This procedure is also implemented without any need for experimental 
testing. ^ 
[0034] Characteristics and advantages of the invention will now be described with reference to an embodiment of 
the invention, illustrated indicatively and by no means exclusively in the accompanying drawings, where: 

30 fig. 1 shows a concentrated-parameter physical model of a tyre used in a method for determining the road handling 

of a tyre of a wheel for a vehicle, constructed according to the invention; 
fig. 2 shows a finite-element model of a tyre used in the method according to the invention; 
figs. 3-6 are schematic representations of test procedures with a non-rolling tyre to which the concentrated-pa- 
rameter physical model of fig. 1 is subjected; 

35 figs. 7-15 are graphs showing the results of the tests illustrated in the figs. 3-6, obtained from the finite-element 

model of fig. 2 which describes a selected real tyre; 

figs. 16-21 depict modes of vibration of the finite-element model describing the selected real tyre; 

fig. 22 is a flow diagram of a procedure for determining the steady state drift curves of the concentrated-parameter 

physical model of fig. 1 ; 

40 fig. 23 is a schematic representation of a brush model associated with the concentrated-parameter physical model 

of fig. 1; 

fig. 24 depicts a contact pressure distribution determined using the finite-element model describing the selected 
real tyre; 

fig. 25 shows details of the brush model of fig. 1 ; 
45 figs. 26 and 27 are graphs illustrating the results obtained with the brush model of fig. 1 ; 

fig. 28 shows schematic representations of the brush model of fig. 1; 

fig. 29 depicts another contact pressure distribution determined using the finite-element model describing the se- 
lected real tyre; 

figs. 30 and 31 are drift curves obtained from the concent rated -parameter physical model describing the selected 
so real tyre; 

figs. 32 and 33 are schematic representations of drift transient state test procedures to which the concentrated- 
parameter physical model of fig. 1 is subjected; 

figs. 34, 35 and 36 are further schematic representations of the brush model of fig. 1 ; 

figs. 37-52 illustrate the drift transient state test results obtained from the concentrated-parameter physical model 
ss describing the selected real tyre; and 

fig. 53 depicts a distribution of forces acting on a beam of the tyre brush model. 

[0035] Illustrated in fig. 1 is a concentrated-parameter physical model, indicated generically with the numeral 1, 
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reproducing the drift behaviour of a tyre of a wheel, made from selected mixes of rubber and reinforcing materials. 
[0O36] The physical model 1 comprises a rigid ring 2 which represents a tread band provided with inserts, a belting 
structure and corresponding carcass portion of the tyre, and a rigid disk 3 representing a hub of the wheel and beading 
of the tyre. The model 1 also comprises principal springsA 5 and 6 and principal dampers 7, 8 and 9 which connect 
the rigid ring 2 to the hub 3 and represent sidewalls of the tyre and air under pressure inside the tyre. The model also 
comprises supplementary springs 10, 11 and 12 and supplementary dampers 13, 14 and 15 which represent phenom- 
ena of deformation of the belt through the effect of a specified vertical load. 

[0037] Associated with the physical model 1 is a brush model 20 which simulates physical phenomena present in 
an area of contact between tyre and road. The brush model 20 comprises a rigid plate 21 under which a system 
representing the tread is applied. The system is preferably bidimensional and comprises numerous parallel, deformable 
beams 22 orientated longitudinally, the ends of which are hinged to the plate, and numerous microinserts, or rows of 
springs. 23 arranged in parallel. In this particular case, there are three deformable beams 22 whilst there are 5 rows 
of microinserts 23 associated with each beam. The bottom ends of the microinserts 23 of the brush model interact with 
a road or the ground 24. The model reproduces the local deformations occurring inside the area of contact and repre- 
sents the uniformly distributed lateral and torsional stiffnesses of the portion of tread in the area of contact. 
[0038] The rigid ring 2 has an equivalent roll radius r [m], mass M c [kg] and diametral moment of inertia J c [kg*m 2 ]. 
The rigid disk 3 has mass M m [kg] and diametral moment of inertia J m [kg*m 2 ]. 

[0039] The principal springs 4 ; 5 and 6 have structural stiffnesses K c , respectively comprising lateral stiffness 
[N/m] between hub and belt, camber torsional stiffness K c9x [Nm/rad] between hub and belt and yawing torsional stiff- 
ness K cHz [Nm/rad] between hub and belt. 

[0040] The principal dampers 7, 8 and 9 have structural dampings R c , respectively comprising lateral damping 
[Ns/m] between hub and belt, camber torsional damping [Nms/rad] between hub and belt and yawing torsional 
damping R cHz [Nms/rad] between hub and belt. 

[0O41] The supplementary springs 10, 11 and 12 have residual stiffnesses K r , respectively comprising residual lateral 
stiffness [N/m], residual camber torsional stiffness [Nm/rad] and residual yawing torsional stiffness [Nm/ 
rad]. 

[0042] The supplementary dampers 1 3, 1 4 and 1 5 have residual dampings Rp respectively comprising residual lateral 
damping R^ [Ns/m], residual camber torsional damping R^ [Nms/rad] and residual yawing torsional damping R^ 
[Nms/rad] 

[0043] The residual stiffnesses and dampings permit allowance to be made for the variation of local stiffness due to 
deflection of the tyre. The lateral and residual yawing stiffnesses and connect the bottom end of the rigid ring 
to the plate, as also does the residual camber stiffness K^. In some cases, the camber deformation of the plate p b 
(absolute rolling rotation) is not taken into consideration, so that connecting the second end of the spring 11, which 
represents the residual camber stiffness, directly to the plate is tantamount to connecting it to the ground. In these 
cases, the stiffness already incorporates the effect due to camber deformability of the brush model. 
[0044] The equivalent system has stiffness per unit length c py and the contact area has a dynamic length 2a and 
dynamic width 2b. 

[0045] With the equivalent system, allowance may be made both for the deformability of the inserts in the tread and 
for the different speeds between a point of the insert in contact with the road (assuming there is adhesion, this point 
has a lateral velocity y' s = 0) and the corresponding point on the belt. Three factors play a decisive role: the coefficient 
of friction at the interface between wheel and road, the normal pressure distribution and stiffness of the inserts in the 
tread. 

[0046] Shown in fig. 1 is an absolute trio of reference axes O-X-Y-2 having versors i, j, k, where the origin O coincides 
with the centre of the hub with the tyre non-deformed, the X axis lies in the plane of the-hub and is of longitudinal 
direction, the Y axis is perpendicular to the X axis and the Z axis is vertical. 
[0047] The degrees of freedom of the physical model 1 are: 

absolute lateral displacement y m of the hub, absolute yaw rotation o m of the hub and absolute rolling rotation p m 
of the hub, 

relative lateral displacement y c of the belt with respect to the hub, relative yaw rotation a c of the belt with respect 
to the hub and relative rolling rotation p c of the belt with respect to the hub, 

absolute lateral displacement y b of the plate, absolute yaw rotation o b of the plate and absolute rolling rotation p b 
of the plate. 

[0048] A further degree of freedom is: 

absolute lateral displacement y s of the bottom ends of the microinserts. 

This degree of freedom has the objective of reproducing the sideward forces created under the contact and which 
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are linked to the relative displacements between the top and bottom ends of the microinserts. In the case of perfect 
. adhesion of the microinserts, with the tyre not drifting, y s = 0. 

[0049] Motion of the physical model is described by assuming small displacements and small rotations of the hub. 

5 [0050] The tyre is described using a Finite-Element Model (F.E.M.) 30 depicted in fig. 2. The finite-element model 
30 comprises first elements (bricks or shells or multilayer composites) with a selected number of nodes, having ap- 
propriately selected form functions, preferably of the first or second order and, even more preferably, linear, and second 
elements suitable for describing the reinforcing materials. Each first element has a first stiffness matrix which is deter- 
mined using a selected characterization of the mixes and a second supplementary stiffness matrix which is determined 

to using a selected characterization of the reinforcing materials. 

[0051] More than ten different types of mixes are usually to be found inside a tyre. Their elastic properties are fed 
into the finite-element model after selected static and dynamic tests are conducted on specimens of the mixes. 
[0052] The static tests consist of tensile, compression and shear tests in which the test conditions, forces and elon- 
gations applied are established in relation to the properties of the mixes measured previously (hardness, etc.). 

15 [0053] The Mooney-Rivlin law of hyperelasticity was taken as the constitutive law. This law describes the specific 
deformation energy in relation to the derivatives of the displacements (deformations), separating the form variation 
energy from the volume variation energy (deviatoric and hydrostatic part of the stress tensor). 

[0054] The coefficients of the constitutive law are calculated in such a way as to minimize the difference between 
the experimental and the calculated deformation energy. 
20 [0055] The dynamic tests are conducted by applying first a static predeformation to the specimens and then an 
oscillating load with frequency in the range from 0.1 to 100 Hz. The dynamic modulus of the mixture is thus detected 
as the complex ratio between stress and deformation. As the frequency is changed, the modulus and relative phase 
between stress and deformation are measured. 

[0056] The reinforcing materials used are fabrics and of metallic type. In the tyres for automobiles, metallic cords 
25 are used only for the belts and fabric cords for the carcass and the outermost belt (zero degrees) located just under 
the tread band. The metallic cord is subjected to pulling until it breaks and to compression in order to obtain by exper- 
imental means the characteristic of the cord from critical load to breaking. This characteristic is implemented in the 
finite-element model. The fabric cord is also subjected to pulling until it breaks. 

[0057] The second elements that describe the reinforcing materials are defined within the first elements (bricks, for 
30 instance) of the finite-element model 30, by assigning the geometrical disposition of the fabrics, the orientation of the 
cord, the spacing between the single cords (thickness) and the experimentally obtained traction and compression, 
characteristic of the cord. These characteristics, in relation to the dimensions of the brick element taken, are resumed 
in a supplementary stiffness matrix which is overlaid on the mixture stiffness matrix, enabling extraction of the cord 
tensions. 

35 [0058] A frequency domain analysis is performed, wherein a linearized response is determined to a harmonic exci- 
tation based on the single degrees of freedom of the physical model. The response is obtained by resolving a matrix 
system, complete with mass, damping and stiffness matrices. Linearization of the matrices is performed at the end of 
the preliminary static determination stages so that non-linear behaviour of the actual tyre is implicitly taken into account. 
[0059] By defining the isotropic linear viscoelasticity, the damping and stiffness matrices in relation to frequency may 

40 be built. It follows that the relation between stresses and deformations is considerably influenced by the elastic behav- 
iour (a higher modulus of elasticity corresponds in particular to higher frequencies) and the damping of the mix. 
[0060] In the case of fig. 2, the finite-element model comprises: 

17,500 elements (16,000 defined by the user + 1 ,500 generated autonomously, needed for definition of the con- 
45 straints of the contact between the tyre and the rim it is fitted on and the road); 

approx. 36,000 nodes (1 9,000 defined by the user + approx. 1 7,000 generated autonomously, needed for resolution 
of the hydrostatic part of the stress tensor); 

approx. 74,000 degrees of freedom (19,000 nodes x 3 Iranslational degrees of freedom for each node + 17,000 
degrees of freedom associated with the elements for resolution of the hydrostatic part of the stress tensor). 

so 

[0061] The equations of motion of the concentrated-parameter physical model 1 of the non-rolling tyre are obtained 
using the Lagrange method. 

[0062] The unknown parameters of these equations are identified by comparing the vibrational dynamic response 
determined using the finite-element model with that obtained from the concentrated -parameter physical model. 
55 [0063] The above-described degrees of freedom (independent variables or generalized coordinates) of the physical 
model are contained in a vector x organized as follows: 
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[0064] The kinetic energy, expressed through the independent variables of the model, is as follows: 



(1.1) 
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(1.2) 



[0065] The potential energy, expressed through functions of the independent variables of the model, is as follows: 
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(1.3) 



where 
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25 



y p is the absolute lateral displacement of the bottom point of the belt; 
s p is the absolute yaw rotation of the belt; 
r p is the absolute camber rotation of the belt. 

The relations linking the physical variables to the independent variables are as follows: 
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p m c 



r p = Pm + Pc 



35 [0066] When these relations are inserted in the potential energy equation, the potential energy in relation to the 
generalized coordinates is obtained. The dissipation energy D is similar in form to the potential energy V 
[0067] If Lagrange's theorem is applied to the kinetic energy expression (1 .2), the inertia terms [M]*x are found, 
where the general mass matrix [M\ is: 
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(1.4) 



[0068] The matrix includes the concentrated parameters of the model M m , M c J m and J c . 

[0069] On derivation of the expression of the potential energy V with respect to the vector of the independent vari- 
55 ables, a general stiffness matrix [K\ is obtained (1 .5): 
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[0070] The structural and residual stiffnesses of the model are included in this matrix. 

[0071] Finally, on deriving the expression of the dissipation energy D with respect to the derivative vector before the 
independent variables, a general damping matrix [R] is obtained (1.6): 
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[0072] The structural and residual dampings of the model are included in this matrix. 

[0073] The vector of the forces F contains the forces and torques acting on the physical model in the different situ- 
ations the model is found in, as will be illustrated later. 
[0074] In short, the equations of motion are obtained in matrix form: 



[M] /x + [R]x + [K]x = F 



(1.7) 



45 [0075] In this equation, the contributions linked respectively to the degrees of freedom of the hub, the belt and the 
tread may be discerned: 
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(1.8) 
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[0076] To identify the concentrated parameters of the physical model, a series of virtual dynamic tests, simulated on 
the finite-element-model 30, is defined, in which the tyre is variously excited and under different conditions so as to 
highlight the contribution of each of the terms that has to be identified. 
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[0077] For each test, the finite-element model is used to determine the frequency response of particular quantities 
measured at precise points of the tyre. 

[0078] To obtain the frequency responses of the concentrated-parameter physical model, the equations of motion 
for the selected tests are defined. 

[0079] For each test and for each frequency response of the quantities taken into consideration (displacements of 
the points considered, constraining reactions, etc.) : the difference is computed between the results obtained from the 
finite-element model and those obtained using the concentrated-parameter physical model. This difference is consid- 
ered an error and is defined as 



[0080] An objective function is then also defined, consisting of the weighted sum of the differences found on each 
channel (each channel corresponds to a quantity considered) between the responses of the concentrated-parameter 
physical model and those of the finite-element model. The objective function may be expressed as follows: 



[0081] Finally each stiffness and damping parameter of the physical model is attributed that value by means of which 
the objective function is minimized, where the vector of the weightings P enables the various channels to be given 
different weightings. 

[0082] To identify the values of the three stiffnesses and the three dampings present between the rigid ring and the 
hub, two virtual tests are conducted on the tyre blown up and not pressed to the ground. Also obtained at this stage 
of the identification are the values for the mass and diametral moment of inertia of the belt that most closely adapt 
behaviour of the model to that of the actual tyre. In particular, the value of the mass of the belt M c is determined in the 
first test and the value of the diametral moment of inertia of the belt J c in the second test. 

[0083] The mass of the hub M m , on the other hand, is determined by imposing the conservation of the total mass of 
the tyre: M tot = M c + M m . A first sharing of the masses is then made, with a view to defining the first-attempt values of 
the mass of the belt and its diametral momeht of inertia. These values will then be modified in the identification stage 
in order to optimize dynamic behaviour of the model. 

[0084] A first test (test A) on the finite-element model consists in exciting the blown-up tyre not pressed to the ground 
by imposing a translation in the Y direction on the hub. The lateral displacement at particular points of the belt of the 
tyre and the force set up between hub and belt are measured. A schematic representation of the first test is shown in 
fig. 3 depicting belt 32, hub 33 in transverse section and a connection spring 34, which represents the lateral stiffness 
and damping between belt and hub. Also shown are an arrow 31 representing the displacement imposed on the hub 
in the Y direction, and a side view of the belt 32, in which the cardinal points N, E, S and W are shown. As the belt is 
symmetrical, the lateral displacement of one cardinal point only of the belt 32 is considered in this test. 
[0085] A first comparison test (A) is carried out on the concentrated-parameter tyre model, as will be illustrated later. 
When the frequency responses of the quantities measured in these first two tests (lateral displacement of a cardinal 
point of the belt and force between belt and hub) are compared and the errors minimized, the following concentrated 
parameters of the model 1 may then be determined: 

the lateral stiffness between hub and belt K cy , 
the lateral damping between hub and belt R cy> 
the mass of the belt M c . 

[0086] A second test (test B) on the finite-element model consists in exciting the biown-up tyre not pressed to the 
ground by imposing a camber rotation on the hub. The lateral displacement at particular points of the belt (camber 
rotation) and the torque transmitted between hub and belt are measured. A schematic representation of the second 
test is shown in fig. 4 depicting the belt 32, the hub 33 in transverse section and a connection spring 35, which sche- 
matically represents the camber torsional stiffness and damping between belt and hub, respectively K cex and R^. 
Also shown are an arrow 41 representing the rotation imposed on the hub about the X axis, and a side view of the belt 
32, in which the cardinal points N, E, S and W are indicated. As the belt is symmetrical, the lateral displacement of 
one cardinal point only of the belt 32, selected indifferently between points N and S, is considered in this test. 



err={v- v FEM ] 



-(1.9) 



N channels 




(1.10) 
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[0087] A second comparison test (B) is carried out on the concentrated-parameter tyre model, as will be illustrated 
later. Again in this case, when the frequency responses of the quantities measured in these second two tests (lateral 
displacement of a cardinal point of the belt and torque between belt and hub) are compared and the errors minimized, 
the following concentrated parameters of the model 1 may then be determined: 

5 

the diametral moment of inertia of the belt J c , 

the camber torsional stiffness between hub and belt K^, 

the camber torsional damping between hub and belt R cex , 

the yawing torsional stiffness between hub and belt K C0Z (by symmetry with K^) 
10 - the yawing torsional damping between hub and belt R c6z (by symmetry with R cex ). 

[0088] The connection between hub and belt is made in like manner for both the camber rotations and the yaw 
rotations and it is possible to consider that the behaviour of the belt not pressed to the ground is axial-symmetrical. 
[0089] A third test (test C) on the finite-element model consists in exciting the tyre blown up, pressed to the ground 

15 and bereft of its tread (rasped) at least in the area of contact by the application to the hub of a force in the Y direction. 
The lateral displacement of the hub and particular points of the belt are measured. A schematic representation of the 
third test is shown in fig. 5 depicting the belt 32, the hub 33 in transverse section, a connection spring 36 (in addition 
to the one described above, indicated with numeral 35) which schematically represents the residual, lateral and camber 
stiffnesses and dampings between belt and ground 24. Also shown are an arrow 51 representing the force imparted, 

20 an arrow 52 representing the lateral displacement of the hub and a front view of the belt 32, in which the cardinal points 
N, E, S and W are indicated In this test, the lateral displacement of the N and E cardinal points of the belt 32 is 
considered. 

[0090] A third comparison test (C) is carried out on the concentrated-parameter tyre model, as will be illustrated later. 
When the frequency responses of the quantities measured in these third two tests (lateral displacement of the hub and 
25 of the N and E cardinal points of the bclU arc compared and the errors minimized, the following concentrated parameters 
of the model 1 may then bo determined: 

the residual lateral stiffness K^. 
the residual lateral damping R^, 
30 - the residual camber stiffness K mx , 
the residual camber damping R rt)x . 

[0091] A fourth test (test D) on the finite-element model consists in exciting the tyre blown up, pressed to the ground 
and bereft of its tread (rasped) at least in the area of contact by the application to the hub of a yawing torque about 

35 the Z axis. The yaw rotation of the hub and the displacement at particular points of the belt are measured. A schematic 
representation of the fourth test is shown in fig. 6 depicting the belt 32, the hub 33 in transverse section, a connection 
spring 37 (in addition to the one described above, indicated with numeral 35) which schematically represents the 
residual yawing stiffnesses and dampings between belt and ground 24. Also shown are an arrow 61 representing the 
yawing torque imparted, an arrow 62 representing the rotation of the hub and a front view of the belt 32, in which the 

40 cardinal points N, E, S and W are indicated. In this test, the lateral displacement of the E and W cardinal points of the 
belt 32 is considered. 

[0092] A fourth comparison test (D) is carried out on the concentrated-parameter tyre model, as will be illustrated 
later. When the frequency responses of the quantities measured in these fourth two tests (yaw rotation of the hub and 
displacement of the.E and W cardinal points of the belt) are compared and the errors minimized, the following concen- 
ts trated parameters of the model 1 may then be determined: 

the residual yawing stiffness K^, 
the residual yawing damping R^. 

50 [0093] During the third and fourth tests, in addition to the identification of the residual quantities, the stiffness and 
damping values already identified in the first and second tests are also modified in order to better describe the transfer 
functions given by the finite-element model. The stiffness and damping values between hub and belt are therefore 
modified slightly from those identified earlier. 

[0094] To obtain the frequency responses of the physical model 1 needed for the identification of its concentrated 
55 parameters, the equations of motion relative to the four tests described above are obtained. 

[0095] The tests A and B are carried out with the tyre blown up and not pressed to the ground, imposing displacements 
on the hub. 

[0096] The free degrees of freedom (generalized coordinates) of the physical model are: 
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whereas the constrained degrees of freedom (constrained coordinates) are: 



^ A ' B = {y m Pm°m} T (1-12) 

[0097] As stated above, only a translation in the Y direction is imposed on the hub in test A and therefore p m = a m 
= 0, whilst only a camber rotation is imposed on the hub in test B and therefore y m = a m = 0. Furthermore the tyre is 
raised off the ground and the residual stiffnesses and dampings relative to the tread cancel each other out. 
[0098] Under these conditions, from the general matrices of mass, stiffness and damping reported above for the- 
concentrated-parameter physical model, a mass matrix [M\ A - B } a stiffness matrix [K\ A ~ B and a damping matrix [R\ A ' B 
are obtained. 

[0099] When a partition is made of the matrices in relation to the free and constrained coordinates, equations of 
motion in scalar form are obtained from (1.8): 



0-13) 



{M m + M c )y m + M e y e =F y 



where F ym , C pm and C am 



represent the forces and torques imposed on the hub. 
[01 00] These equations govern the motion of the system in the simulated tests A and B. 
[0101] For the displacements imparted to the hub, the equations are: 



<T_ = 



- J <P,n 



(1.14) 



[01 02] The frequency responses must now be determined of the quantities measured in the tests A and B to be able 
to compare the frequency responses given by the finite-element model with the frequency responses of the concen- 
trated-parameter model. 

[0103] In the test A, with the physical model, the frequency response of the lateral displacement of the belt is recon- 
structed, taking into -consideration the first equation of the system (1.14). The displacement imparted on the hub is: 



(1.15) 



the complete differential equation to be resolved is as follows: 



KPZ+J^y/ +K„y c A = M e .Cl\A.e 



(1.16) 



[0104] The solution of this equation, having a sinusoidal force agent equivalent to y m , is of the type: 



y c A = B.e int (1 . 17) 



[0105] By means of substitutions and simplifications, the frequency response is obtained for the degree of freedom 
y c upon variation of Q, that is to say the frequency response of the lateral displacement of the belt with respect to the hub. 



12 



EP 0 949 496 A1 



A M c& A M 

y c = 2 



2 m 

e /iif (1.18) 



[0106] As the finite-element model provides the absolute displacement of the points of the belt, the frequency re- 
sponse of the absolute lateral displacement of the belt may be determined from the absolute lateral displacement: 

• 2 
10 AAA. M Q. 1 , . /O/ M -|Q\ 

y c _ass =y c + y m = i t - » — a — — + 1 i Aa (119) 

[0107] Upon computing the difference between the lateral displacement of the physical model and that provided by 
is the finite-element model, an error is obtained that is a function of M c , and R cy . In the identification stage, this error 
is minimized, that is to say, values are chosen for M c K cy and R cy that make the error minimum. 
[0108] In the test A, the frequency response of the force transmitted between hub and belt provided by the finite- 
element model is also compared with that provided by the concentrated-parameter physicaf rnodel. This force is equal 
to: 

20 

F^S^+ISLRJaJ™ (1.20) 

[0109] Again with this quantity, when the difference is computed between the frequency response of the force de- 
25 termined with the finite-element model and that obtained using the concentrated-parameter physical model, a second 
error is obtained that is also function of M c and R^. The above-stated quantities are obtained upon minimization 
of this error. 

[0110] In the test B, with the physical model, the frequency response of the rotation of the belt is reconstructed and 
the same method of procedure as the test A is then followed. The second equation of the system (1 . 1 4) is considered, 
30 starting from a camber rotation imparted to the hub of the type: 

Pm e =C.e /£2 ' 0-21) 
35 the frequency response is obtained of the absolute camber rotation of the belt: 

B = a B + 1]C .e ,m (1.22) 

repass r c K m L 2 . ^ „ is \ 

(-J e .Q +i.a.R e0x + K eix ) 

40 

and, from this, the frequency response of the rotation of the lateral displacement of the S point of the belt is 
determined: 

45 v B =P B r (1.23) 

[0111] Also determined is the frequency response of the torque transmitted between hub and belt: 

C hub ^ lt B = (K oex + i.CX.R c6x ).D.e int (1.24) 

[0112] Upon computing the difference between the displacement y c a ss_south 8 and the torque C hub . beU of the physical 
model and those given by the finite-element model, an error is obtained that is a function of the quantities J c , K C0X and 
55 | n t h e identification stage, this error is minimized and values are identified for the above-mentioned quantities. 

[0113] In the tests C and D, with the tyre blown up, pressed to the ground and bereft of tread, the free degrees of 
freedom and the constrained degrees of freedom of the physical model are determined in the two test situations and 
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an axial force is imposed on the hub in the Y direction, in the test C, and a yawing torque about the Z axis, in the test 
D. When these conditions are put into the equations of motion (1.7) and (1.8), the equations of motion for the free 
degrees of freedom of the physical model are obtained and the complex vector of the free degrees of freedom is 
determined. This vector contains both the modulus and the phase of the free degrees of freedom in question. 
[0114] In the test C, from this vector, the frequency responses of the lateral displacement of the hub, of the N point 
and of the E point of the belt are determined. In the identification stage, the difference is computed between the fre- 
quency response obtained from the finite-element model and that obtained from the concentrated-parameter physical 
model and a second error is obtained which is a function of K^, R^, and R^. On minimizing this error, the above- 
mentioned quantities are determined. 

[0115] In the test D, from the above-mentioned vector, the frequency responses are determined of the yaw rotation 
of the hub and of the displacement of the E point of the belt in relation to the quantities and R^. These quantities 
are determined by minimizing the error resulting as the difference between the frequency response provided by the 
finite-element model and that provided by the concentrated-parameter physical model. 

[0116] Once the concentrated parameters of the physical model have been identified, it is possible to determine 
is some quantities that allow an evaluation, or an initial approximation at least, to be made of the behaviour of the tyre. 
These quantities are the total drift stiffness of the tyre K d , comprised in turn by the structural stiffness K c the tread 
stiffness K b and the total camber stiffness Ky of the said tyre. 

[0117] The stiffness K d is an important parameter for the purposes of definition of the tyre in that it provides an 
indication of the effect that the design parameters have on the drift behaviour. 
20 [0118] To obtain K d the free and constrained degrees of freedom are partitioned: 



25 
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35 



40 



45 
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[0119] A further partition is made of the vector x, into internal (of the belt) and external (of the tread) degrees of 
freedom: 

?# = <* e ?/J 7 = &o Pc a c» T (1 26) 

[0120] The stiffness matrix (1.5) is then modified following these partitions and the equations of motion can be written 
as follows for the free degrees of liberty: 

Ifr]*. +[*.]*, -£,=s < L27 > 

[0121] Remembering that no external forces act on the degrees of internal freedom, the second equation of the 
system (1.27) may be resolved in relation to the external degrees of freedom, the result inserted in the first equation 
and the latter made explicit in relation to the external coordinates x e alone. At this point a stiffness matrix is defined 
for the tyre structure, limited to the external degrees of freedom: 

I &=-WKk (l ' 28) 

[01 22] The resolving equation is: 

[£*c]x< =Z e (1.29) 



wherein: 



DKicnnrin- ^ctd na^a-aaAA 1 i 
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(1.30) 



is a matrix (2,2). The first element of the first row of this matrix represents the total stiffness of the tyre K d : 



K d = 



(1.31) 



15 



[0123] With the stiffnesses of the tyre with concentrated parameters known from the identification procedure, it is 
possible to calculate K d . 

[0124] The total camber stiffness Ky of the tyre provides an indication of the tyre's ability to exit from a longitudinal 
track namely parallel to the direction of movement of the tyre, made on the road when the tyre is moving with rectilinear 
motion The method according to the invention enables calculation of Ky without the need to make and test a tyre 
prototype, as is generally done. 

[0125] To determine K^, a virtual test practically equal to the experimental one, is carried out consisting in imposing 
a selected camber angle on the hub, with a null drift angle, and in measuring the sideward force created on the hub. 
[0126] The free and constrained degrees of freedom are identified: 



20 



(1.32) 



25 



ij={y m Pm °m Pb> 



(1.33) 



[01 27] By imposing a camber angle and leaving all the other constrained degrees of freedom unaltered, the sideward 
force is determined by resolving the matrix equation: 



30 



[K].x = f 



(1.34) 



35 



where the matrix [K\ is the (1.5) and F a vector containing the forces acting on the tyre. Figuring among these 
are the force and torque due to the deformations undergone by the inserts of the tread in the assumption that there is 
perfect adhesion. These two contributions are: 



40 



F yb = 2c p' a ~y b 
as 2 



(1.35) 



45 



[0128] As the previous two terms are functions of two of the nine degrees of freedom of the concentrated-parameter 
model (y b and o b ), they must be expressed in function of these degrees of freedom and thus be brought to the left- 
hand side of the equals sign in the matrix equation reported above (1 .34), so that they too contribute to determining 
the stiffness matrix of the system. 

[0129] The partitioned matrices K^, K, v v , K vl ? and K w ? are known: 
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(1.36) 
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(1.37) 



(1.38) 
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(1.39) 

30 [01 30] From these, the vector of the forces acting on the constrained degrees of freedom may be calculated using 
the equation: 



(1.40) 



[0131] Remembering that: 



40 



pj=-Ki-[^r.[^r-HKj 



(1.41) 
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is a matrix (4,4), the first element of the second row, which is the total camber stiffness of the tyre, is determined: 

(1.42) 



K y = #(1,2) 



[01 32] To demonstrate the validity of the method according to the invention, the results obtained are brought into the 
procedure for the identification of the concentrated parameters of the physical model of a tyre of the 55 range (H/C 
section ratio of 0.55) manufactured by the Applicant. In the graphs of figs. 7-15, the frequency response obtained with 
the concentrated-parameter physical model is represented with a continuous line, while the frequency response ob- 
tained with the finite-element model during the tests A, B, C and D is represented with asterisks. 
[01 33] Fig. 7 is a graph of the frequency response of the displacement of a point of the belt obtained during the test 
A from the physical model and from the finite-element model, whereas fig. 8 is a graph of the frequency response of 
the force created between the hub and the belt, again in the test A. 

[0134] The values of the concentrated parameters identified with the test A are reported in the following Table I. 
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Parameter 


Symbol 


Value identified 


Linear stiffness between hub* and belt 9 


K cy 


5 e 5 [N/m] 


Linear damping between hub and belt 


R cy 


116.906 [Ns/m] 


Belt mass 


Mo 


7.792 [kg] 



[01 35] It may bexibserved from the graphs of figs. 7 and 8 that the frequency response has a single resonance peak. 
Corresponding to this peak is a mode of vibration of the tyre illustrated in fig. 16. This mode of vibration maintains the 
tread band substantially rigid and may therefore be described by the concentrated-parameter model, that simulates 
the sidewall and the belt with a rigid ring. The parameters identified are therefore valid in a set of frequencies ranging 
from 0 to 100 Hz, since modes of vibration that considerably deform both sidewall and belt appear with frequencies 
that are higher than this. 

[0136] Fig. 9 is a graph of the frequency response of the displacement of a point of the belt obtained during the test 
B from the physical model and from the finite-element model, whereas fig. 10 is a graph of the frequency response of 
the torque created between the hub and the belt, again in the test B. 

[01 37] The values of the concentrated parameters identified with the test B are reported in. the following Table II . 

Table II 



Parameter 


Symbol 


Value identified 


Torsional stiffness between hub and belt 


K cex /K c6z 


3.3 e 4 [N/m] 


Torsional damping between hub and belt 


R ce>/ R cez 


8.217 [Ns/m] 


Diametral moment of inertia of belt 


Jc 


0.373 [kg+m 2 ] 



[0138] As in the previous test, again in this one there is a single resonance peak corresponding to the modal defor- 
mation of the tyre shown in fig. 17. 

[01 39] In the test C, consideration also needs to be given to the normal load bearing on the tyre. Three loads corre- 
sponding to three standard working conditions were applied to the tyre: a reduced load of between 2,500 and 3,000 
N, an intermediate load of between 3,500 and 4,800 N and a high load of between 5,100 and 6,500 N. 
[0140] Fig. 11 is a graph of the frequency response of the lateral displacement of the hub obtained from the physical 
model and from the finite-element model during the test C, with a normal load of 2,914 N bearing on the tyre, whereas 
figs. 12 and 13 are graphs of the frequency response of the lateral displacements of the E and N points of the belt, 
again in the test C. 

[0141] The values of the concentrated parameters identified in the test C are reported in the following Table III. 

Table III 



Parameter 


Symbol 


Value identified 


Residual linear stiffness 




665964 [N/m] 


Residual camber stiffness 


K t8x 


11461 [Nm/rad] 


Residual linear damping 


R ry 


2042.345 [Ns/m] 


Residual camber damping 


R r6x 


0.358 [Nms/rad] 



[0142] Unlike in the two earlier tests, there are two resonance peaks in this third test corresponding to two modes 
of vibration of the tyre depicted in figs. 18 and 19. 

[0143] Again in this case, the modes of vibration maintain the sidewall and belt complex substantially rigid and may 
therefore be described accurately by the concentrated-parameter model, in other words the concentrated-parameter 
model is valid in. the range of frequencies between 0 and 100 Hz. 

[01 44] The frequency responses obtained with the loads of 4,611 N and 6,302 N are not shown herein. The results, 
however, are similar to those shown for the 2,914 N load. 

[0145] Fig. 14 is a graph of the frequency response of the yaw rotation of the hub obtained from the physical model 
and from the finite-element model during the test D, with a normal load of 2,914 N bearing on the tyre, whereas fig. 15 
is a graph of the frequency response of the lateral displacement of the E point of the belt, again in the test D. 
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[0146] The values of the concentrated parameters identified with the test D are reported in the following Table IV. 



Table IV 



Parameter 


Symbol 


Value identified 


Residual yaw stiffness 


K r6z 


14933 [Nm/rad] 


Residual yaw damping 




14.318 [Nms/rad] 



[01 47] Again in the test D, there are two resonance peaks which correspond to the two modes of vibration of the tyre 

10 shown in figs. 20 and 21 . 

[01 48] In this case, whereas the first mode of vibration (fig. 20) is comparable to a rigid mode, it is much more difficult 
to describe the second one (fig. 21 ) as a rigid mode. It seems apparent that the belt as a whole becomes deformed in 
this second mode of vibrating. The reconstructed frequency response does not therefore accurately reproduce the 
second resonance peak. The concentrated-parameter model in the test D is therefore only valid in the range of fre- 

75 quencies from 0 to 70 Hz. 

[0149] Again in the test D, the frequency responses obtained with the loads of 4,611 N and 6,302 N are not shown 
herein The results, however, are similar to those shown for the 2,914 N load. 

[0150] To evaluate a tyre in relation to its road handling, verification needs to be made that.the total drift stiffness K d 
of the tyre and the total camber stiffness KL are within the following value ranges: 

20 * 

K d =500 - 2,000 [N/g] 

2S K y = 40 - 3,500 [N/g] 

K c = 8,000 - 30,000 [N/g] 

30 

K b = 150 -400 [N/g] 

where g = degree. 

[01 51] The method according to the invention enables an evaluation to be made of the steady state behaviour of the 

35 drifting tyre. The bidirectional brush model shown in fig. 1 and illustrated above is used. 

[0152] Under drift conditions, the microinserts of the brush model become deformed and a sideward force and a 
moment of torque act on the beam that they are connected to. These forces and moments result in a deformation of 
the beam that affects the configuration of the microinserts. By means of successive iterations, the deformation effec- 
tively assumed by the beams in relation to the drift angle imposed on the hub is determined. At this point, the total 

40 sideward force is determined, together with the total self-aligning torque acting on the plate that the beams are con- 
nected to. Under the total sideward force and self-aligning torque, the plate snakes by an angle that depends on the 
overlying structure, i.e. on the springs connecting it to the hub. On account of this snaking, the deformation of the , 
microinserts is altered. By performing new iterations, the sideward force and self-aligning torque that the tyre summons 
up in reaction to the drift angle imposed are determined. 

45 [0153] The procedure is illustrated in the flow diagram of fig. 22. 

[0154] A drift angle a is applied to the physical model 1 , starting from a condition in which the beams 22 are in a 
non-deformed configuration and the brush model has a null snaking o b (block 45). The sideward forces acting on the 
beams through the effect of the drift angle a and depending on the difference a - o b and on the deformation of the 
beams are determined (blocks 46a, 46b and 46c). The deformation of the beams is determined (blocks 47a, 47b and 

50 Ale). A check is made to see if the deformation is the same as that determined in the previous step (blocks 48a, 48b 
and 48c). The procedure for determining deformation of the beams is repeated until deformation is verified to be equal 
to that found in the previous step. At this point, the snaking of the plate, i.e. of the brush model, is determined (block 
49). A check is made to see if the snaking is equal to that calculated in the previous step (block 50). The procedure 
for determining snaking of the plate is repeated until snaking is verified to be equal to that found in the previous step. 

55 At this point, the sideward force and the self-aligning torque acting on the hub due to the drift angle imposed are 
determined (block 51 ). The procedure is repeated for the different values of the drift angle a to produce drift, force and 
self-aligning torque curves that enable steady state drift behaviour of the tyre to be evaluated. 
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[0155] The beams of the brush model do not all become deformed in the same way since the pressure distribution 
acting on each beam is different, as also are the sideward forces. In practice, the sidewalls of the tyre (shoulders) have 
greater stiffness than in the central band of the tread. 

[0156] To identify the flexural stiffness of each beam, a distributed parameters model is used based on a model 

5 depicted in fig. 23 wherein the equivalent beams and springs (microinserts) of the brush model are indicated with the 
same numerals as in fig. 1. A linear static analysis, corresponding to a sideways traction and made by applying a 
sideward force to the contact surface of the road, is then conducted on the finite-element model. The accordingly 
stressed tyre becomes deformed and the lateral displacements are determined at circumferential sections thereof, on 
a level with the external band. These circumferential sections are divided into three groups corresponding to the central 

10 part of the area of contact and to the two parts at either side. The mean lateral displacement of each group is calculated 
over the full length of the contact area. On determining the difference between the generic deformation mean (mean 
lateral displacements over the full length of the contact area) provided by the finite-element model and the lateral 
deformation of the rigid ring ol the concentrated-parameter physical model, also subjected to the same sideward force 
as applied to the finite-element model, a "difference 0 deformation is obtained that must be offset by the equivalent 

15 beam of flexural stiffness EJ (N*m 2 ). 

[0157] Knowing how the sideward force is divided over the three parts of the contact area from the linear static 
analysis of sideways traction on the finite-element model, a flexural stiffness EJ is determined for each beam. 
[0158] In fig. 23 it may be seen that the sections closest to the sidewall that the traction is exerted on become more 
deformed than those further away (dashed lines). If a load of F/3 were to be applied to each beam, there would be a 

20 monotone pattern of E J when moving from one side of the tyre to the other, but this would be contrary to what experience 
shows. Therefore, using the finite-element model, a horizontal division of the load between the various sections is also 
made. 

[0159] The results of the identification procedure are reported in the following Table V. 



Table V 



30 



EJ [N^m 2 ] 


1 st load 


2 nd load 


3 rd load 


1 st beam 


50 


120 


250 


2 nd beam 


45 


106 


210 


3 rd beam 


50 


120 


250 



[0160] To describe as accurately as possible what happens in the area of contact between tyre and road, the actual 
pressure distribution in this area must be determined. 

[0161] Using the finite-element model, a static pressure distribution is determined on a non-rolling tyre with no drift. 
The pressure distribution determined under a 55 range tyre, quoted earlier, is illustrated in fig. 24, with a normal acting 
load of 2,914 N. Similar distributions were determined with normal loads of 4,611 N and 6,302 N. 
[0162] It may be seen that two pressure peaks are presented in the transverse direction and that these peaks are 
shifted outwardly with respect to the centre of the area of contact. Furthermore the distribution in the longitudinal 
direction is symmetrical. Finally the ratio of the maximum value to the minimum value of the pressure along the trans- 
verse direction increases as the normal load increases. 

[0163] To obtain the drift curves, the pressure values in those points where the microinserts of the brush model are 
present must be known beforehand. The contact was represented in the cortcentrated-parameter physical model by 
means of a regular grid of 200 x 1 5 elements (200 elements longitudinally and 3x5 elements transversally). Accordingly, 
the pressure at the nodes of this grid must be known, i.e. at 3,000 points. The finite-element model provides the pressure 
at a much lower number of points arranged in an irregular grid. A procedure for the interpolation of the finite-element 
model data is used to move from the grid of points provided by the finite-element model to that required by the con- 
centrated-parameter model. In the case in hand, the procedure adopted is that of the "inverse distances". This inter- 
polation procedure requires the following inputs: 

the coordinates of those points at which the pressure is known; 
the value of the pressure at these points; 

the coordinates of the points for which the interpolated pressure value is required. 
[0164] The procedure provides the pressure value for the points required. 

[0165] Downstream of the identification stage, it is required that any negative pressure values be equal to zero. 
Accordingly the actual shape of the area of contact can be reproduced extremely accurately as each microinsert of 
the brush model has an own length and a known pressure distribution. More specifically, the length of the contact of 
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each microinsert changes upon variation of the transversal position considered and each microinsert has an own 
pressure value that depends on its position within the contact. Therefore the curve of the contacts of each microinsert, 
which identifies the position of the bottom ends of the microinserts when the tyre is drifting, varies from one microinsert 
to the next. When the wheel is set rolling and drifting, and a beam and the five rows of microinserts under it are 
considered, it may be seen that the microinserts undergo a progressive, linearly increasing deformation. The range of 
deformations is therefore triangular, in the passage from entrance to exit of the footprint area, assuming there is no 
slipping (fig. 25). Fig. 26 illustrates the deformation undergone by the beam due to the effect of the sideward forces 
transmitted by the microinserts and the deformations undergone by the microinserts. The deformations of the micro- 
inserts are given by the relative displacements between the upper ends (attached to the beam) and the lower ends 
(located on the line of contacts in the event of adhesion). Pattern of the pressure in each microinsert is illustrated in 
fig- 27. 

[0166] In order to determine the lateral and torsional stiffness per unit of length of the microinsert, respectively c p 
and k.tor tne total stiffnesses of the whole tread are taken and then the values found are divided by the total length of 
the microinserts obtained as the sum of those of the single microinserts: 
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2.a.c„ 



2.a.k„ 



(1.43) 



where 2a is the length of the area of contact, l tot is the sum of the lengths of the single microinserts, c py is the 
25 stiffness per unit of length of the contact and k tor is the torsional stiffness per unit of length of the contact. 
[0167] These values are used to obtain the drift curves. 

[0168] To determine the torsional contribution of the microinserts in the yaw test with a non-rolling tyre, the flexure 
of the inserts present in the footprint area of the concentrated-parameter model is identified, the sum is obtained of 
the contributions of each and the value obtained subtracted from the torque per foot provided by the finite-element 
30 model. The shape of the area of contact is provided by the finite-element model. In this way, the purely twisting torque 
arising following rotation about the Z axis of the tyre pressed to the ground is obtained. 

[0169] The sideward force produced following deformation of each micrpinsert and contributing to the flexure, as- 
suming perfect adhesion, is equal to the lateral stiffness of the microinsert multiplied by the relative displacement in 
the plane of the contact of the top end with respect to the bottom end of the microinsert. Assuming a rotation about 
35 the centre of contact and perfect adhesion, the deformation of a generic microinsert is equal to d (x, y) * a', where a' 
is the rotation about the Z axis effectively undergone by the tread and d is the distance of the microinsert taken from 
the centre of the contact. 

[0170] The lateral stiffness of the single microinsert is determined starting from the stiffness per unit of length of the 
contact: 

40 

c p = ^&.dx.dy (1.44) 

45 where a is the half-length of the contact, b is the half-width of the contact and dx and dy are respectively the 

longitudinal and transversal dimensions of each microinsert. 
[0171] The bending torque due to each microinsert is calculated as follows: 

so M f (x, y) = F(x, y). d{x, y) (1 .45) 

where F is the sideward force due to the flexure of the microinsert and x and y are the coordinates as measured 
from the centre of the contact. 

[0172] The total bending torque created under the area of contact for the rotation a' is calculated as follows: 

55 
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* / jr y x y 

s [01 73] By subtracting the total bending torque from the torque per foot provided by the finite-element model, the pure 
twisting torque is obtained: 

M^C^-M, 0-47) 

10 

[0174] If the twisting torque is divided by the rotation a' and by the length of the area of contact 2a, the torsional 
stiffness per unit of length of the contact is obtained: 

M 

is K = tors (1.48) 

[0175] With the method according to the invention making allowance for the deformability of the structure of the 
concentrated-parameter model, of the local deformations of the contact through the equivalent beams, of the lateral 
20 and torsional stiffnesses of the microinserts, the drift curves for the range 55 tyre were determined for three different 
loads. 

[0176] To obtain calculated drift curves coinciding with the experimental ones, the shape of the area of contact is 
considered. 

[0177] Bearing in mind the local deformation of the belt in correspondence with the area of contact, a pear-shaped 
25 area of contact is identified and a brush model with equivalent beams of different lengths is used, as depicted in fig. 
28. The length of the beams is obtained by performing a sensitivity analysis, assuming that the central beam is equal 
in length to the statistically measured contact area length and that the variation of length of the external and internal 
beams is of equal modulus, but opposite sign. The sensitivity analysis was conducted on the range 55 tyre. 
[0178] It was found that the shape of the area of contact affects the transversal pressure distribution: the pressure 
30 peakcorresppnding to the outside of the curve is no longer equal to the other one (the external pressure peak increases, 
whereas the internal one decreases as the drift angle increases). It is considered that the pressure distribution in the 
transverse direction is always less symmetrical with a higher pressure peak in correspondence with the longest equiv- 
alent beam and a lower pressure peak in correspondence with the shortest equivalent beam. 

[0179] For an indication of the pattern of the distribution of pressure from the finite-element model, the sideward 
3S forces arising under the contact are calculated for different drift angles through the concentrated-parameter physical 
model as described earlier. This force distribution is applied to the nodes of the finite-element model (non-rolling) in 
contact with the ground and an extremely realistic indication of the pressure distribution in the area of contact obtained. 
From the pressure distributions obtained, a strong symmetry in the transverse direction results. Illustrated in fig. 29 is 
the pattern of pressure distribution (N/mm 2 ) for a vertical load of 2,914 N whenever a sideward force and a self -aligning 
torque corresponding to a test with drift of 6° are applied to the finite-element model. On application of these new data 
to the concentrated-parameter physical model, calculated drift, sideward force and self-aligning torque curves are 
obtained that are practically coincident with those obtained by experimental means. Figs. 30 and 31 illustrate respec- 
tively the pattern of the sideward force (N) and the self-aligning torque (Nm) in relation to the drift angle a (°) for a 
vertical load of 2 : 914 N. The calculated values are marked with asterisks (*), whereas the experimental values are 
45 marked with circles (o). Similar curves were obtained for normal loads of 4,611 N and 6,302 N. 

[0180] To determine the dynamic behaviour of the drifting tyre during the transient state, the tyre is made roll by 
imposing laws of motion on the hub that vary with time and are suitable for numerically reproducing selected experi- 
mental tests commonly carried out on tyres. 

[0181] Two experimental tests are conducted in the laboratory for evaluating the dynamic behaviour of a drifting tyre: 



50 



55 



a first test called pendulum test, consisting in simultaneously imposing a lateral displacement and a steering angle 
on the hub of the tyre; 

- a second test called drift test with yaw pattern, consisting in directly imposing a drift angle equal to the yaw angle 
imparted to the hub. 

[0182] In both of the experimental tests, a wheel-road that simulates the ground is used. The tyre is blown up and 
pressed to the wheel-road which is rotating at constant angular speed. The axis of the tyre is made oscillate, and a 
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drift angle that varies with time induced upon it. The oscillations are imposed on the wheel with the two test arrange- 
ments described above. 

[0183] Fig. 32 illustrates a cinematic diagram of a test machine with the pendulum arrangement. By adopting a 
reference system integral with the hub of the wheel, the longitudinal and transverse components are defined of the 
velocity of the centre of the area of footprint determined by the rotation velocity V imparted to the tyre by the wheel-road. 

V T = &J + V.sin{<j) 

V L =F.cos(cr) (L49) 

from which the drift angle a is obtained for small values of the angle of pendulum (steering) a: 

as aJ+_V^o = a+ a / (150) 

[0184] If a motion that is, for instance, sinusoidal is imposed on the pendulum: 

G = a 0 .cos(Q.f) (1-51) 

the drift angle becomes: 

a = a 0 .cos(n.f) - — ^./.s/n(Q-0 = a 0 .(cos(Q.f) + <p) (1-52) 

[0185] Therefore, the drift angle a is equal to the steering angle a imposed on the hub, with a phase difference of 
angle <p. If a is constant, then a is always = o. 

[0186] Fig. 33 illustrates a cinematic diagram of a test machine with the yaw arrangement. The longitudinal and 
transverse components of the velocity of the centre of the hub are defined as follows: 



V L =V.cos(a) - (1.53) 

V T = V.sin(cr) 



[0187] The resultant drift angle for low snaking angles imposed on the hub is: 

a = -£H&Z!=o (154) 
, V.cos(a) 

[0188] Again in this case the steering angle imposed on the hub is equal to the drift angle. 

[0189] To determine the equations of motion of the tyre under dynamic conditions, suitable for reproducing the two 
experimental tests described above, the absolute, right-handed trio of reference axes shown in fig. 1 is adopted. The 
horizontal X axis, orthogonal to the axis of rotation of the hub, forms a null steering angle am. In the numerical simu- 
lations it is assumed that it is the hub that moves while the wheel-road remains motionless. The microinserts of the 
brush model enter the area of contact and their bottom ends remain attached to the ground until slipping occurs. 
[0190] The independent variables of the physical model 1 with nine degrees of freedom are again those indicated 
above (1.1). The equations of motion governing the motion of the physical model are again those in matrix form reported 
above (1 .7). The matrix of mass and that of stiffness are those reported above at (1 .4) and (1 .5). The damping matrix 
is that given above at (1 .6), wherein the terms (5,3), (5,6) and (6,5) are equal to J ycoy (moment of inertia and angular 
speed about the Y axis) to take gyroscopic effects into account. 
[0191] The vector of the forces comprises: 

the force F by and the torque M bz l transmitted by the microinserts to the structure of the tyre and a function of time; 
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the resulting sideward force F by is positive if orientated similarly to the Y axis and the self-aligning torque M b2 is 
positive if orientated similarly to the Z axis; 

the force F my and the torques M mx and M m2 transmitted by the test machine to the hub; 
- the torque M bx constraining the plate not to rotate about the X axis because such a degree of freedom is not 
5 included since a monodimensional brush model is used. 

f = {? my M mx M mx 0 0 0 F by M bx M bz ] T d .55) 

io [01 92] For the case under examination, in the concentrated-parameter physical model in the transient state, the yaw 
- test is simulated by stating: 

ymSSPa Z? (1.56) 

[0193] Then the three free coordinates of the hub are constrained in the model whereas the roll of the brush model 
is not taken into consideration since the contact is considered to be monodimensional. 
[0194] A partition is made between the generalized x, and constrained x^degrees of freedom: 



20 



?= {*?s v V= Pc °c y b *bHy m Q m * m o b }} r (1.57) 

and the matrices of mass, stiffness and damping are accordingly rearranged by performing a division into four 
25 subrnatrices so as to obtain two matrix equations from the equations of motion: 



| [M u ]x, + [M lv ]x v + [R u ]x, + [R lv ]x v + [K u ]x, + [K lv ]x v = F t 



(1.58) 



[01 95] In the first equation, F, represents the vector containing the external active forces acting oh the actual degrees 
of freedom; in this case : the only external forces acting on the x, are the sideward force F by and the self -aligning torque 
35 M bz created under the contact and acting on the degrees of freedom of the plate. These forces depend on the arrange- 
ment of deformation of the inserts and on whether or not there is local slipping; they are therefore, generally speaking, 
non-linear functions of the degrees of freedom. In the second equation, F^ represents the vector of the generalized 
reactions applied to the constrained degrees of freedom. 

[0196] Remembering that the vector x v and its derivatives with respect to time are vectors of known functions, the 
40 first equation may be rewritten as follows: 

[M u ]x, + [/?„]*, + [K u ]x, =F,- [M h ]x v - [R h ]x v - [K,„]x„ = F (1 .59) 

45 

wherein the terms F_ are all known because they are the sum of the effective external forces and of the equivalent 
forces due to the motion imparted to the constraints. 

[0197] The equation (1.59) therefore represents a system of "n u equations, one for each unknown x,; by solving this 

equation, the motion X| of the model can be determined. 
so [0198] Once the equations (1.59) have been integrated and the values obtained for x, and their derivatives, the 

constraining reactions F v may be obtained from the equations for the constrained degrees of freedom. 

[0199] To solve the non-linear equations (1.59) with explicit numerical methods, the procedure is to invert the matrix 

of mass [M f j\ and, as the matrix is singular and cannot therefore be inverted, the equation (1 .59) is rewritten as a first 

order system. A further partition is then made, followed by a change of variables. 
55 [0200] The selected partition of x, (1 .57) is: 
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[0201] On rearranging the already rearranged matrices of mass, stiffness and damping as described above, and 
performing a division into four submatrices, the equation (1 .59) becomes: 



\[M ce ]k c +fc c ]x c+ [R cb ]x b+ [K cc ]x c+ [K cb ]x b =F C 



(1.61) 



where F C/ is a vector of three elements containing the generalized forces acting directly on the degrees of freedom 
of the belt and F b is a vector of two elements containing the generalized forces acting directly on the degrees of freedom 
of the brush model (F by and M b2 ). 

[0202] Finally, to obtain a first order system, an auxiliary identity is added to the system (1 .61 ): 



[MccK + K]i £ + K]x b + [K ee ]x c + [K cb ]x b = F c 
[M cc ]x c = [M a ]x c 



(1.62) 



and the following change of variables is made: 



[0203] If the matrices [B[ and [C] are defined: 



(1.63) 



lB} = 



[c]= 



K] [0] [R cb ] 
[0] [0] fcj 
. [0] [M cc ] [ojj 

' M fc] M 

M fc] K] 

-K] [o] [o] 



(1.64) 



the system (1.62) is synthetically expressed as: 



[B}z + [C}z = F 



(1-65) 



where F is a vector of eight elements comprised as follows: 



E = {£/ tj o o o} r 



(1.66) 



[0204] The matrix [B\ is now invertible. These dynamic equations are numerically integrated a Runge-Kutta step-by- 
step method of the 3 rd order. 
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[0205] In the equations (1 .65), the matrices [B\ and [Q are defined whilst the vector of the forces still needs to be 
identified, in particular F by and M bz created under the rolling tyre in transient state. 

[0206] As already stated earlier, to describe the behaviour of the rolling tyre in the transient state, a component of 
displacement V*t in the X direction is imposed on the centre of the hub. In the drift test with yaw arrangement, small 

5 variations of the degrees of freedom are imposed on the physical model. Under these conditions, the deformation 
modalities of the microinserts inside the area of contact differ from those under steady state conditions. An insert enters 
the area of contact in a generic position that depends on the motion of the wheel and its deformation in the contact 
area changes with time in a manner dictated by the pattern of the model's degrees of freedom. The forces arising 
during the transient state depend on the trajectory of the top and bottom ends of each microinsert and are, as stated, 

10 influenced by the motion of the wheel. The deformation of a generic i-th microinsert is defined as follows: 



Y U<** l) - Y Lp( x i J) (1 ' 67) 

15 where Yj. P is understood to be the absolute lateral displacement of the top end of the i-th microinsert connected 

to the plate and Y;_ T is understood to be the absolute lateral displacement of the bottom end interacting with the ground. 
[0207] A procedure enabling an instant-by-instant evaluation of the displacements of the top and bottom ends of all 
the microinserts in the contact is now illustrated with reference to figs. 34 and 35. The plate 21 of fig. 1 , connecting the 
belt to the microinserts in a generic configuration, is illustrated in schematic form, seen from above, in fig. 34. Having 

20 earlier indicated the absolute lateral displacement of the plate with y b and its absolute yaw with o* b , the absolute lateral 
displacement of the top end of the generic microinsert Y,.p is : 



Y j p (x, t)=Y b (t) + (a- sin(o b ) (1 .68) 

where ^ is the abscissa indicating the generic microinsert in the reference integral with the wheel (fig. 34). 
[0208] The position must now be identified of the bottom ends of the microinserts Y hT in order to determine their 
deformation and accordingly obtain the contact forces. 

[0209] A generic microinsert is taken at the generic instant of time at which it enters the contact area: the top and 
bottom ends, seen from above, are at the same point because the microinsert has not yet been deformed. In the 
successive instants, the top end of the generic microinsert has a longitudinal displacement component V * t opposite 
to the feed direction of the tyre while, at the same time, due to the effect of the absolute lateral displacement of the 
plate y b and of its yaw o b , it also has a transverse direction component. Simultaneously, the bottom end of said micro- 
insert, assuming perfect adhesion, remains motionless in the absolute reference system whereas, in the reference 
system integral with the wheel, it possesses a displacement component V * t in the longitudinal direction. Having set 
a distance between two adjacent microinserts of V * dt in correspondence with each step of integration of the dynamic 
equations described above, a single microinsert enters and exits from the contact (fig. 35). Three successive instants 
of time are taken. At the instant t, the microinsert 1 has just entered under the contact and is in the deformed config- 
uration. An instant later, the top end of the microinsert 1 has gone into 1 ' having a vertical coordinate £ = V * dt, whereas 
the bottom end has remained in 1 . Simultaneously, the non-deformed microinsert 0 comes under the contact. The 
position occupied by the bottom end of the microinsert 1" (which is the third microinsert in the local reference) at the 
time t + 2 * dt is the same as occupied by the bottom end of the microinsert V an instant earlier In the same way, the 
position occupied by the bottom end of the microinsert 2", the fourth in the local reference, at the time t + 2 * dt is the 
same as that occupied by the microinsert 2' an instant earlier. The position occupied by the bottom end of the generic 
microinsert at the time t is: 



Y L7 <x if t)= Y LJ (x._^ t-dt) (1.69) 

[0210] The non-linearity of the brush model results from the possibility of whether the microinserts slip or not. Two 
factors have a fundamental role in this: the coefficient of friction at the interface between the wheel and ground and 
the pressure distribution. A side view is shown of the reference model, at the bottom of fig. 36. The pressure distribution 
is assumed to be parabolic in the longitudinal direction, whereas no pattern need be defined for the transverse direction 
since the brush model considered is monodimensional. 

[0211] The normal forces per unit of length that are discharged to the ground are: 



BNSDOCIDkEP 0949496A1 L> 



25 



EP 0 949 496 A1 



where F z is the vertical load applied to the tyre. The maximum deformation possible for the generic microinsert is: 

where u is the coefficient of friction between tyre and ground and c py is the stiffness per unit of length of the brush 
model. When this maximum deformation is brought onto the line x representing the top ends of the microinserts, the 
area within which the bottom ends must fall is determined (seen from above, top part of fig. 36). Knowing the position 
of the top end and the bottom end of the generic microinsert, its deformation can be determined. 
[0212] When the deformation of the microinserts and the stiffness of the tread are known, jt is possible to determine 
the forces created under the contact area in the Y direction. 
[0213] The resulting sideward force F by acting on the plate is: 

n 

^(') = Z(^_r(*i -.0- !!_,(*„ 0-<>.Ax (1.72) 



[0214] If the single sideward forces created under each single microinsert are integrated and multiplied by the re- 
spective arm, the self-aligning torque M bz is found as follows: 

M k (0 = Y J (Y iT (x i ,t)-Y iP (x i ,t).c I ,.Ax.x j (1.73) 

[0215] Some results obtained for a rolling tyre from the procedure described are illustrated below. 
[0216] To simulate the behaviour of the tyre in the transient state, a test was conducted envisaging a step input of 
the steering angle imparted to the hub. This test is particularly important as it serves to evaluate the time taken by the 
tyre to go into a steady state and therefore, in practice, it gives its speed of response. 

[0217] In this test, a selected steering angle is imposed instantaneously on the hub at a time t - 0 and the dynamic 
behaviour of the tyre is obtained from the time history of the six free degrees of freedom of the model, of the sideward 
force and of the self-aligning torque. 

[0218] The test was conducted on a range 55 tyre with an applied vertical load of 2,914 N and a steering angle of 
four degrees imposed on the hub. The feed speed of the tyre was 30 Km/h. 

[0219] Fig. 37 shows the time pattern of the lateral displacement of the belt y c with respect to the hub and of the 
absolute lateral displacement of the plate y b . 

[0220] In figs. 38, 39 and 40, the time pattern is shown respectively of the snaking a c of the belt, of the snaking c b 
of the plate and of the roll p c of the belt. 

[0221] Figs. 41 and 42 illustrate the time pattern respectively of the sideward force F b and of the self-aligning torque 
Mbz- 

[0222] The pattern of the sideward force F by is used to evaluate the length of relaxation which consists of the space 
travelled by the tyre before the sideward force reaches 63.2% of its steady state value. 

[0223] The behaviour of the tyre in drifting in the transient state is reproduced using, as the transversal drift force 
applied to the vehicle F by that given by the following equation: 

%rF by+ F by =F by (a) (1.74) 

where a is the instantaneous drift angle, F by is the steady state force and 5 is the length of relaxation. 
[0224] In order to determine the length of relaxation, the signal of the sideward force F by is removed of its high 



26 



45 



50 



EP 0 949 496-A1 

frequency harmonics. A low-pass filter was used with cutoff frequency 30 Hz (fig. 43). The filtered signai is comparable 
to the response of a 1 st order system receiving a step signal as its input. 
[0225] The general equation of motion of such a system is as follows: 

A.f t + B.F=F extemgl (1.75) 

which, through the Laplace transform, may be written as follows: 

(1.76) 



10 ~ 2 

F 1 B a 



External As+B £ S+ 1 XS+1 
B 

is where t is the time constant of the system and represents the time taken by the system to reach 63.2% of its 

steady state value. The (1 .76) is the system transfer function. If a step of amplitude Xq is provided as input, the system 
response in the frequency domain is as follows: 

F(s)4 0 -|^ (1.77) 
and : on passing to the time domain, the following is obtained: 

25 



F(t) = a.x 0 ^l-e ') (1.78) 

30 [0226] The term a * x 0 is the steady state value of R 

[0227] When the equations (1.74) and (1.75) are compared, the length of relaxation is obtained from the curve of 
the filtered sideward force. 

[0228] From the filtered signal of fig. 43, the steady state value of the sideward force is obtained so as to calculate 
63.2% and. from here, trace back to the time constant x. For an immediate evaluation of the steady state value, the 
35 same filtering operation was also carried out on the self -aligning torque (fig. 44). 

[0229] In the case in question, the relaxation length was found to be 0.23 metres. 

[0230] Figs. 45-47 illustrate the results obtained with the range 55 tyre for three different loads and with different yaw 
angles imposed on the hub with feeding speed of 100 Km/h. 

[0231] Local deformability of the tyre in the transient state is also introduced into the brush model by means of the 
40 equivalent beam described above. After determining the deformation curve, the positions of the top ends of the micro- 
inserts are identified whereas the position of the bottom ends is determined using the procedure described above. On 
the basis of the data obtained, the single elementary forces acting on each section of a beam can be computed (fig. 53). 
[0232] The resulting sideward force and the self -aligning torque are determined by integrating the forces distributed 
along the beam. The constraining reactions F a and F b are determined in order to give the deformation curve. 



a 

F y = \q,*dx (1.79) 

—a 
a 

Mz = \q y *x*<£c (1.80) 
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(1.81) 

3 2 2*a K } 
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b 2 2*a 



(1.82) 



[0233] 



Curvature of the beam is related to the bending torque through the following equation: 



EJy"(x)=+M£x) 



(1.83) 



[0234] This equation is resolved by determining the bending torque in each section of the beam. 
[0235] Figs. 48, 49 and 50 illustrate the pattern respectively of the sideward force, of the self-aligning torque, of the 
length of relaxation at steady state in relation to the steering angle, with a rigid and deformable contact, for the range 
55 tyre to which a vertical load of 2,914 N is applied. 

[0236] Figs. 51 and 52 illustrate respectively the pattern of the length of relaxation and of the time constant in relation 
to velocity for the range 55 tyre to which three vertical loads of 2,914 N, 4,611 N and 6,302 N were applied. 



1 . Method for determining the road handling of a tyre of a wheel for a vehicle, said tyre being comprised by selected 
mixes of rubber and reinforcing materials, said method comprising: 

a) a first description of said tyre by means of a first, concentrated-parameter, physical model (1), said first 
physical model (1) comprising a rigid ring (2) which represents the tread band provided with inserts, a belting 
structure and corresponding carcass portion of said tyre, a disk (3) which represents a hub of said wheel and 
beading of said tyre, principal springs and dampers (4, 5, 6, 7, 8, 9) connecting said rigid ring (2) to said hub 
(3) and representing sidewalls of said tyre and air under pressure inside said tyre, supplementary springs and 
dampers (10, 11, 12, 13, 14, 15) representing deformation phenomena of said belting structure through the 
effect of a specified vertical load and a brush model (20) simulating physical phenomena in an area of contact 
between said tyre and a road, said area of contact having a dynamic length 2a ; 

b) a definition of selected degrees of freedom of said first physical model (1), and 

c) an identification of equations of motion suitable for describing the motion of said first physical model (1 ) 
under selected dynamic conditions, 

characterized in that it comprises 

d) the definition of said concentrated parameters, said concentrated parameters consisting of the mass M c 
and a diametral moment of inertia J c of said rigid ring, the mass M m and a diametral moment of inertia J m of 
said disk, structural stiffnesses K c and structural dampings R c respectively of said principal springs and damp- 
ers, and residual stiffnesses K r and residual dampings R r respectively of said supplementary springs and 
dampers, wherein 



- said structural stiffnesses K c consist of lateral stiffness K cy between said hub and said belt, camber 
torsional stiffness between said hub and said belt and yawing torsional stiffness K cez between said 
hub and said belt, " 

- said structural dampings R c consist of lateral damping R cy between said hub and said belt, camber 
torsional damping R c9x between said hub and said belt and yawing torsional damping R c02 between said 
hub and said belt, 

- said residual stiffnesses K r consist of residual lateral stiffness K^, residual camber torsional stiffness 

and residual yawing torsional stiffness K^, and 

- said residual dampings R r consist of residual lateral damping R^, residual camber torsional damping 
R^ and residual yawing torsional damping R rt2 , 

e) a description of said tyre by means of a second, finite-element model (30) comprising first elements with a 
selected number of nodes, suitable for describing said mixes, and second elements suitable for describing 
said reinforcing materials, each first finite element being associated with a first stiffness matrix which is deter- 
mined by means of a selected characterization of said mixes and each second element being associated with 
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a second supplementary stiffness matrix which is determined by means of a selected characterization of said 
reinforcing materials, 

f) a simulation on said second, finite-element model (30) of a selected series of virtual dynamic tests for exciting 
said second model (30) according to selected procedures and obtaining transfer functions and first frequency 
responses of selected quantities, measured at selected points of said second model (30), 

g) a description of the behaviour of said first physical model (1 ) by means of equations of motion suitable for 
representing the above dynamic tests for obtaining second frequency responses of said selected quantities, 
measured at selected points of said first physical model, 

h) a comparison between said first and said second frequency responses of said selected quantities to deter- 
mine errors thai are a function of said concentrated parameters of said first physical model (1), 

" and 

i) the identification of values for said concentrated parameters that minimize said errors so that said concen- 
trated parameters describe the dynamic behaviour of said tyre, 

j) the determination of selected physical quantities suitable for indicating the drift behaviour of said tyre, and 
k) the evaluation of the drift behaviour of said tyre by means of said physical quantities. 

2. Method according to claim 1 , characterized in that said selected physical quantities are the total drift stiffness K d 
of said tyre, in turn comprising the structural stiffness K c and the tread stiffness K b , andlhe total camber stiffness 

of said tyre. 

3. Method according to claim 1 , characterized in that it also comprises 

c) a definition of said brush model (20), said brush model (20) having a stiffness per unit of length c py and 
comprising at least one rigid plate (21 ), at least one deformable beam (22) having a length equal to the length 
2a of said area of contact and at least one microinsert (23) associated with said beam (22), said microinsert 
(23) consisting of at least one set of springs distributed over the entire length of said beam (22), said springs 
reproducing the uniformly distributed, lateral and torsional stiffness of said area of contact. 

4. Method according to claims 1 and 3, characterized in that said degrees of freedom referred to at previous point b) 
are composed of: 

absolute lateral displacement y m of said hub, absolute yaw rotation o m of said hub and absolute rolling rotation 
p m of said hub, 

relative lateral displacement y c of said belt with respect to said hub, relative yaw rotation o c of said belt with 
respect to said hub and relative rolling rotation p c of said belt with respect to said hub, 

absolute lateral displacement y b of said plate, absolute yaw rotation a b of said plate and absolute rolling rotation 

p b of said plate, 

and 

absolute lateral displacement y s of the bottom ends of said at least one microinsert. 

5. Method according to the claim 1 , characterized in that said selected series of virtual dynamic tests referred to at 
previous point f) comprises a first and a second test with said tyre blown up and not pressed to the ground, said 
first test consisting in imposing a translation in the transverse direction y on the hub and in measuring the lateral 
displacement y c of at least one selected cardinal point of said belt and the force created between said hub and 
said belt in order to identify said mass M c said lateral stiffness K cy> and said lateral damping R^, said second test 
consisting in imposing a camber rotation 6* on said hub and in measuring the lateral displacement of at least one 
selected cardinal point of said belt y c and the torque transmitted between said hub and said belt in order to identify 
said diametral moment of inertia J c , said camber torsional stiffness K cex , said camber torsional damping R c0x> said 
yawing torsional stiffness K c92 and said yawing torsional damping R^. 

6. Method according to claims 1 and 5, characterized in that said selected series of virtual dynamic tests referred to 
at previous point f) also comprises a third and a fourth test with said tyre blown up, pressed to the ground and 
bereft of said tread at least in said area of contact, said third test consisting in applying to said hub a sideward 
force in the transverse direction F y and in measuring the lateral displacement y c of said hub and of at least two 
selected cardinal points of said belt in order to identify said residual lateral stiffness K^, said residual lateral damping 
R^, said camber residual stiffness K^, and said camber residual damping R^, said fourth test consisting in 
applying to said hub a yawing torque C ez and in measuring the yaw rotation of said hub and the lateral displacement 
y c of at least one selected cardinal point of said belt in order to identify said residual yawing stiffness and said 
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residual yawing damping R^z- 

7. Method according to claims 1 and 3, characterized in that it also comprises 

m)an application to said first physical model (1) of a drift angle a, starting from a condition in which said at 
least one beam (22) is in a non-deformed configuration and said brush model (20) has a null snaking a b> 
n) the determination of the sideward force and the self-aligning torque that act on said hub through the effect 
of said drift and which depend on the difference a - a b and on the deformation of said at least one beam (22), 

0) the determination of the deformation curve of said at least one beam (22), 

p) an application of said sideward force and said self-aligning torque to said second, finite-element model (30) 
in order to obtain a pressure distribution on said area of contact and 

q) the determination of the sideward force and the self-aligning torque that act on said hub through the effect 
of said drift a on said first physical model (1 ), that depend on the pressure distribution calculated in the previous 
step p), 

r) a check, by means of said pressure distribution obtained in the previous step p), that said sideward force 

and said self-aligning torque are substantially similar to those calculated in previous step q), 

s) a determination of the sideward force and of the self-aligning torque for said angle of drift, and 

t) repetition of the procedure from step m) to step s) for different values of the drift angle a to obtain drift, force 

and self-alignment torque curves, suitable for indicating the drift behaviour under steady state conditions of 

said lyre, and 

u) the evaluation of the steady state drift behaviour of said tyre. 

8. Method according to claim 1 , characterized in that it also comprises 

1) a simulation of the behaviour of said first physical model (1 ) in the drift transient state by means of equations 
of motion reproducing selected experimental drift tests, and 

ii) the determination, with a selected input of a steering angle imposed on said hub, of the pattern with time 
of the selected free degrees of freedom of said first physical model, of the sideward force and of the self- 
aligning torque in said area of contact in order to determine the length of relaxation of said tyre. 

9. Method according to claim 1 , characterized in that said first elements of said second, finite-element model (30) 
have linear form functions and their stiffness matrix is determined by means of selected static and dynamic tests 
conducted on specimens of said mixes, whereas the stiffness matrix of said second elements is determined by 
means of selected static tests on specimens of said reinforcing materials. 

10. Tyre for a wheel of a vehicle, said tyre being made from selected mixes of rubber and reinforcing materials and 
comprising a carcass, a belting structure, a tread band provided with inserts, shoulders, sidewalls, beads provided 
with bead wires and bead fillings, said tyre being representable by means of a first, concentrated-parameter, phys- 
ical model (1) and a brush model (20) with a road, characterized in that said concentrated parameters comprise 
structural stiffnesses K c consisting of lateral stiffness K cy , camber torsional stiffness and yawing torsional 
stiffness K c02 , structural dampings R c consisting of lateral damping R cy , camber torsional damping R c6x and yawing 
torsional damping R^, residual stiffnesses K r consisting of residual lateral stiffness K^, residual camber torsional 
stiffness and residual yawing torsional stiffness K^ z , and residual dampings R r consisting of residual lateral 
damping R^, residual camber torsional damping R^ and residual yawing torsional damping R^, 

said tyre also being representable by means of a second, finite-element model (30) comprising first elements with 
a selected number of nodes, suitable for describing said mixes, and second elements suitable for describing said 
reinforcing materials, said concentrated parameters being identified by means of a selected series of dynamic 
tests on said second, finite-element model and represented by equations of motion applied to said first physical 
model, said tyre having construction characteristics substantially equivalent to said concentrated parameters which 
describe the dynamic behaviour of said tyre and enabling the determination of selected physical quantities suitable 
for indicating the drift behaviour of said tyre for evaluation of said tyre in relation to its road handling. 

11. Tyre according to claim 10, characterized in that said selected physical quantities are the total drift stiffness K d of 
said tyre, in turn comprising the structural stiffness K c and the tread stiffness K b , and the total camber stiffness Ky 
of said tyre. 

12. Tyre according to claim 11 , characterized in that the total drift stiffness K d and the total camber stiffness K y are 
within the following value ranges: 
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K d = 500 - 2,000 [N/g] 

K = 40 - 3,500 [N/g] 
5 7 

where g = degree. 

1 3. Tyre according to claim 11, characterized in that the structural stiffness K c and the tread stiffness K b are within the 
io following value" ranges: 

K c = 8,000 - 30,000 [N/g] 

75 K b = 150 -400 [N/g] 

where g = degree. 
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Test no. 1: rim lateral displacement 
TO i ! ! ! ! ! r 




O 10 20 30 <Q SO 60 70 BO 90 TOO 

freq. [Hz] 




freq. [H2] 



F\ C 7 



Test no. 1 : reaction on rim 
1 OOOO i ! * ! : 




freq. (Hz) 
f \ Cr. 8 



35 



EP 0 949 496 A1 
Test no. 2: lateral displacement North 
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Test no. 3, load no. 1: rim lateral displacement 
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Test no. 3, load no. 1: belt lateral displacement East 
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Test no. 3, load no. 1 : belt lateral displacement 
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Test no. 4, load do. I : belt lateral displacement East 
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Pressure distribution 
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Sideward force with time 
3SOOi 1 — , — : 




FIG. UA 




RM.Qnrv-tm- <-pd noAQ^Ai i -> 



50 



EP 0 949 496 A1 



Sideward force with time 
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Steady-staie self-aligning torque with variation of ihe vertical load 
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Steady-state self-aligning torque with rigid and deformable contact 
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